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ABSTRACT 

Observations of the southern Cepheid £ Car to yield the mean angular diameter and 
angular pulsation amplitude have been made with the Sydney University Stellar In- 
terferometer (SUSI) at a wavelength of 696 nm. The resulting mean limb-darkened an- 
gular diameter is 2.990±0.017mas (i.e. ±0.6 per cent) with a maximum-to-minimum 
amplitude of 0.560±0.018mas corresponding to 18.7±0.6per cent in the mean stellar 
diameter. Careful attention has been paid to uncertainties, including those in mea- 
surements, in the adopted calibrator angular diameters, in the projected values of 
visibility squared at zero baseline, and to systematic effects. No evidence was found 
for a circumstellar envelope at 696 nm. The interferometric results have been com- 
bined with radial displacements of the stellar atmosphere derived from selected radial 
velocity data taken from the literature to determine the distance and mean diameter 
of I Car. The distance is determined to be 525±26pc and the mean radius 169±8i?Q. 
Comparison with published values for the distance and mean radius show excellent 
agreement, particularly when a common scaling factor from observed radial velocity 
to pulsation velocity of the stellar atmosphere (the p-factor) is used. 

Key words: Cepheids: general — stars: individual (£ Car, HR 3884) — stars: dis- 
tances — techniques: interferometric 



1 INTRODUCTION 

The accurate determination of the zero point of the Cepheid 
Period-Luminosity relation is an important step in the re- 
finement of the extragalactic distance scale. The combi- 
nation of interfe rometric and spectroscopic data, as dis- 
cussed by iDavid (1999), is an approach to this task that 
is now feasible. The first demonstration of the combination 
of interferometric and spectroscopic data to determine both 
th e mean diameter and distance of a Cepheid was made 
bv lLane et al] \200dt ) based on interferometric observations 
with the Palomar Testbed Interfer o meter (PTI) of ^ Gem. 
iLane. Creech-Eakman fc Nordgr"enl (|2002l ) subsequently ap- 
plied the technique to rj Aql and, with addi tional data, to 
C Gem , again based on PTI observations. iKervella et al.l 
l|2004al ) have reported the measurement of the angular pul- 
sations of seven Galactic Cepheids made with the European 
Southern Observatory's Very Large Telescope Interferom- 
eter (VLTI) and, in combination with spectroscopic data, 
have determined the distances to the Cepheids. 

One of the prime programmes for which the Sydney 
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University Stellar Interferometer (SUSI) (jPavis et al.lll999l ') 
was developed was the measurement of the mean diameters 
and angular pulsation amplitudes of Cepheids but, in its 
initial configuration, it lacked the required sensitivity. The 
recent commissioning of a red beam-combination system in 
SUSI (|Davis et all l2007al ) has resulted in a significant in- 
crease in sensitivity. In this paper we report SUSI measure- 
ments of the mean angular diameter and angular pulsation 
amplitude of £ Car, one of th e seven Cepheids measured with 
the VLTI bv lKerveUa et all |2004a ) . A preliminary analysis 
was presented a t a European Sout hern Observatory sym- 
posium in 2005 (|Davis et al.ll2007bl ') and IJacobI l|2008l ) has 
discussed the SUSI Cepheid programme including an anal- 
ysis of a subset of the £ Car observations reported here. 

Spectroscopic radial velocity measurements from the lit- 
erature have been combined to produce a radial velocity 
versus pulsation phase curve and, after scaling with a fixed 
p- factor as discussed in Section [6. 3 1 this has been integrated 
to give the radial displacement of the Cepheid surface as 
a function of pulsation phase. The radial displacement and 
limb-darkened angular diameter data have been combined 
to determine the distance and mean radius of £ Car. 

The resulting distance and mean radius are compared 
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with values in the hterature determined by the same tech- 
nique and by the infrared surface brightness method and 
exceUent agreement is found, particularly when all values 
are scaled to a common p- factor. 



2 THE PULSATION PHASE 

Although it might seem premature to discuss the pulsation 
phase prior to discussing the observations it is necessary to 
ensure that both interferometric and spectroscopic phases 
are synchronised since the epochs of the interferometric and 
spectroscopic observations do not overlap. A problem in 
the case of I Car is the fact th at there i s evid ence that 
its period has undergone changes. ISzabadoi ([l98^) reported 
that the period of I Car changed from 35.5318 days prior to 
JD2440000 to 35.5513 days. Shobbrook (1992) refined the 
Szabados period , with additional photome try in 1990, to 
35.5443 days and berdnikov fc Turneil (|2004l l. based on pho- 
tomet ry in 2002, give a value of 35.5572 days. iTavlor et al.l 
have listed values back to 1901 and found in their 
own spectroscopic study of i Car from early 1991 to late 
1996 that an increase in the period was evident although 
the uncertainties were large. 

The approach that we have adopted is as follows. Since 
the largest body of radial velocity data available to us 
is the combined Mount Stromlo Observatory (MSO) and 
Mount John University Observa tory (MJUO) observations 
analysed by iTavlor et a l.' (^19971) a nd th eir phases are ef- 
fectively based on the fshobbroqj (|l992t ) ephemeris, albeit 
with small adjustments to bring the early MJU O data into 
line w ith the MSO data, we have adopted the IShobbrookl 
l|l992h ephemeris for the analysis of both interferometric 
and additional spectroscopic data. The adopted period is 
35. 5443±0. 0006 days with 2447880.81±0.10 the zero point 
of maximum light in heliocentric Julian Date. The phases of 
individual data points have been computed from the Julian 
Date of observation for both interferometric and spectro- 
scopic data. 

Because of the variations in period it is found that there 
are small phase shifts between data sets taken at different 
epochs. The change in phase across any given data set as 
a result of adopti ng the Shobbrook p e riod r ather than for 
example, that of .Berdnikov fc Turned (|2004h . is negligible, 
and a small phase shift of the whole data set to bring it 
into alignment with the ITavlor et al.l (|l997l ) data is justi- 
fied. Details of the alignment of the different sets of data 
will be discussed when combining data and, in particular, 
spectroscopic data in Section |6] 



3 THE INTERFEROMETRIC OBSERVATIONS 

Measurements of the squared frin ge visibility V'^ were 
made with SUSI l|Davis et al.l Il999l ) using the red beam- 
combinati on system, which employs the fringe-scanning 
technique l|Davis et al.ll2007ah . This beam-combination sys- 
tem uses matched filters with a central wavelength of 700 nm 
and spectral bandwidth of 80 nm. The observations to de- 
termine the angular diameter of I Car were made with a 
baseline of 40 m with additional measurements at 5 m to en- 
able the zero baseline value of to be checked. 



Table 1. Calibrators used for the observations of I Car. mggs 
is the estimated magnitude at 695 nm, 9uD the adopted uniform 
disk angular diameter at 695 nm and f2 the angular distance of 
the calibrator from I Car. 
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3.1 Calibrators 

Calibrators were selected as close in the sky as possible to 
I Car with the additional requirement of being minimally 
resolved. The limiting sensitivity of SUSI at 700 nm is ~+5 
which limited the choice of calibrators and compromise was 
necessary. The calibrators used are listed in Table [1] with 
their spectral types, adopted uniform-disk angular diame- 
ters, and angular distances from i Car. Common calibrators 
were used throughout to eliminate the potential influence 
of calibrator diameters on the pulsation curve. The effective 
wavelength of observations of all three calibrators has been 
estimated to be 6 95.0±2.0 nm f ollowin g a similar analysis to 
that described bv lDavis et al] l|2007ah . 

The uniform-disk angular diameters have been 
determined from measurements made with the 
Narrabri Stellar Intensity Interferometer (NSII) 
(|Hanburv Brown. Davis fc AUenlll974l ) and with the Mark 
III O ptical Interferometer (Mark III) (|Mozurkewich et al.l 
120031 ). In the case of /3 Car the value measured with the 
NSII has been adopted after correction from the limb- 
darkened angular diameter to the uniform-disk angular 
diameter using the appropriate corr ection factor for 695 nm 
interpolated from the tabulation of iDavis. Tango fc BoothI 
(2000). In the absence of measured angular diameters for 
t and s Car, limb-darkened angular diameters have been 
determined by interpolation in a plot of NSII and Mark 
III limb-darkened angular diameters for unreddened visual 
magnitude Vo = 0.0 versus {B-V)o for (B-V)o < 0.6. The 
interpolated limb-darkened angular diameters have been 
corrected to the Vb magnitudes of the stars and corrected for 
limb-darkening to give the uniform-disk angular diameters 
listed in Table [1] again using appr opriate corre c tion f actors 
interpolated from the tabulation of lDavis et al.i (|200G| ). The 
uncertainties in the uniform-disk angular diameters have 
been estimated from the scatter in the plot of limb-darkened 
angular diameters versus {B-V)i3- 

3.2 The Observations 

Observations of calibrators and I Car were alternated in 
each observing session so that every £ Car observation was 
bracketed by observations of calibrators. Each observation 
consisted of a set of 1000 scans each 140 ^m long and con- 
sisting of 1024 by 0.2 ms samples. Each scan set was followed 
by photometric and dark scans. One complete observation 
of I Car, including the bracketing calibrators, took a total 
of ~18 minutes. 

Observations were made with a baseline of 40 m on 31 
nights for (. Car between 2 March 2004 and 24 May 2007. 
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Initially all 31 nights were included in the analysis but four 
nights were subsequently rejected when their values were 
found to lie more than four standard deviations from the fit 
to limb-darkened angular diameter versus radial displace- 
ment of the stellar surface (Section [7]) and the analysis was 
then repeated without them. The maximum deviation from 
the fit after rejection of the four points was < 3.1a. Ex- 
amination of the data for three of the four rejected nights 
revealed that they were obtained during a short period in 
late March-early April 2004 when there was significant leak- 
age of the metrology laser light into the signal beams. On 
the fourth rejected night, the data were poor, showing large 
scatter with only two out of the six points lying within one 
standard deviation of the mean for the night. 

The projected baseline (i.e. the effective baseline of an 
observation) was less than 40 m due to the southerly decli- 
nations of £ Car and its calibrators. 

Observations were also made with a 5 m baseline, close 
to the phases of maximum and minimum angular diameter, 
on 13 March and 5 May 2007. These will be discussed in 
Section 15.11 



4 ANALYSIS OF INTERFEROMETRY 

The initial analysis of the fringe scans for £ Car and the 
calibrators was car ried out in the SUSI software "pipeline" 
l|Davis et al.ll2007al ') which outputs the raw and seeing cor- 
rected squared visibility amplitudes [V'^), the Julian Date, 
projected baseline, hour angle, fluxes etc. for each set of 
scans. For the early observations up to 7 April 2004 full sets 
of photometry files were not recorded. In these cases pho- 
tometry files were duplicated to enable the scans to be pro- 
cessed in the pipeline and it is noted that analysis has shown 
that this procedure has a negligible effect on the resulting 
values of . For each night the output file from the pipeline 
was imported into an Excel spreadsheet for examination of 
the data and for further analysis. 

For each set of calibrator scans a Transfer Function (T) 
was calculated. T is given by 



Table 2. Calibrated values of V^. The format for the date (Uni- 
versal Time) is ddmmyy, JD is the mean Julian Date of the ob- 
servations minus 2450000, b is the mean projected baseline, and 
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5 THE ANGULAR DIAMETER 

The uniform-disk angular diameter (^ud) for each night was 
determined by fitting the equation 



y cxp 



(1) 



where is the observed, seeing corrected value of and 
Kxp is the expected value of calculated from the uniform- 
disk angular diameter, effective wavelength and projected 
baseline. 

The values of Vq^s for ^ Car were then multiplied by the 
weighted mean value of T for the two calibrators bracketing 
the £ Car observation to give the 'true' calibrated value of 
for £ Car. That is 



,(£Car)= ""^^ 



Vl,{£ Car) 



T 



(2) 



The mean value of Krue(^ Car) for each night is listed 
in Table [2] with the date (Universal Time), the mean Ju- 
lian Date ( JD) , the mean projected baseline, the calibrators 
used, and the number of scan sets for £ Car. The uncertainty 
in V^y^e{£ Car) takes into account the uncertainties in the 
uniform-disk angular diameters of the calibrators although 
these are generally negligible compared to the uncertainties 
in the values of for both calibrators and £ Car. 



2Ji(x) 



(3) 



to the individual values of Vtruc(^ Car) assuming that 
Vtruc(^ Car) at zero baseline was unity. In equation Q Ji{x) 
is a Bessel function and x = nb 6vu /A where b is the pro- 
jected baseline and A is the wavelength of the observation. 

The effective wavelength for £ Car will be a function of 
pulsation phase as the star changes colour during its cycle. 
Calculations of effective wavel ength for supergian ts follow- 
ing the procedure described bv lDavis et al.l ()2007ah give the 
mean effective wavelength for £ Car as 696.2 nm with a vari- 
ation of ~ ±0.7 nm, or ±0.1%, over a complete pulsation 
cycle. The systematic uncertainty in the calculated effective 
wavelength is conservatively estimated to be ±2.0 nm and 
we therefore adopt a fixed value of 696.2±2.0 nm as the ef- 
fective wavelength for £ Car. Since the expected variation 
is significantly less than the adopted uncertainty and, gen- 
erally, small compared with the observational uncertainties 
in the angular diameters, the adoption of a fixed value is 
justified. 
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5.1 The Zero Baseline Value of Vii^^{i Car) 

It is important to ciieck tlie validity of the assumption 
that Vtruc(^ Car) is unity at zero baseline. If it differs from 
unity as a result of the presence of a companion or scat- 
tered light from extended circumstellar matter, for exam- 
ple, it would translate to errors in the angular diameters. 
A large circumstellar envelope (CSE) has been discovered 
ar ound I Car at A^-band using long-baseline interferometry 
bv lKervella et al.l (|2006l ) and its signature was detected by 
the same authors in their A'-band observations. CSEs have 
als o been detected in th e A-band around Pol aris and 5 Cep 
bvlMerand et all ||2006|') and a round Y Oph bv lMerand et al.l 
l|2007^ . iNardetto et al.l (|2008h confirm a dominant absorp- 
tion component in Ha for i Car, whose velocity is constant 
and near zero in the stellar rest frame, and they attribute 
this to a CSE. These detections of CSEs, particularly around 
I Car, emphasise the importance of checking the validity of 
our assumption. In order to carry out this check, observa- 
tions were made at a baseline of 5 m on 13 March 2007 at 
phase 0.023 and on 5 May 2007 at phase 0.514. The ob- 
servations were therefore close to minimum and maximum 
angular diameter respectively. The observations were made 
after the initial analysis of the 40 m data had been completed 
so values for the uniform-disk angular diameter at the two 
phases could be interpolated from a plot of uniform-disk 
angular diameter versus phase. Given the angular diame- 
ter, and assuming T4rue(^ Car) to be unity at zero baseline, 
the expected value of Krue(^ Car) at the projected base- 
line for each night was calculated. If the assumption that 
^rue(^ Car) is unity at zero baseline is valid the observed 
and expected values will be in agreement. The results are 
listed in Table [3] and show that at minimum angular diam- 
eter the observed and expected values of K?uc(^ Car) dif- 
fer by 0.004±0.013 and at maximum angular diameter by 
0.004±0.012. The assumption that Vt^^^{t Car) is unity at 
zero baseline is therefore justified and there is no evidence 
of a significant effect of a CSE at 696 nm. 



5.2 The Uniform-Disk Angular Diameters 

The uniform-disk angular diameters corresponding to the 
dates of observation, together with the corresp onding pul- 
sation phases determined using the ephemeris of lShobbrookl 
as discussed in Section [2l are listed in Tabled 
Figure[T]shows the data and fitted uniform-disk angular 
diameter curves for near minimum angular diameter (phase 
= 0.047 on 24 May 2007) and near maximum angular diam- 
eter (phase = 0.503 on 6 February 2005). 



5.3 The Uniform-Disk Angular Diameter of the 
Calibrator t Car 

A check was made on the consistency of the transfer function 
T determined for the three primary calibrators, namely /3, 
t and s Car. For each pair of consecutive observations of 
calibrators the transfer functions were plotted against each 
other. For example, where observations of (5 Car and s Car 
bracketed an observation of I Car, the transfer function of 
/3 Car was plotted against the transfer function of s Car. For 
all the 40 m observations this resulted in 67 data points for 
s Car V. P Car, 85 points for l Car v. P Car, and 51 points 



Table 4. The angular diameter of £ Car. JD is the mean Julian 
Date of the observations minus 2450000 (see Table|2j, <^ the mean 
pulsation phase for the observations, SuD the uniform-disk angu- 
lar diameter, Red. the reduced of the fit of equation ((Sjl 
to the values of V^^^^{i Car), P696 the limb-darkening factor, and 
9ld is the limb-darkened angular diameter. Further details are 
given in the text. 
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Figure 1. Examples of the data and fitted uniform-disk angular 
diameter curves for observations near minimum and maximum 
angular diameters. Key: The filled circles are the data for phase 
0.047 on 24 May 2007 and the solid line is the fit to these data; 
the open circles are the data for phase 0.503 on 6 February 2005 
and the dashed line is the fit to these data. Details of the fitting 
procedure are given in the text. 
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Table 3. The results of the 5 m observations near minimum and maximum values of the angular diameter of £ Car. The format for the 
date is ddmmyy, JD is the mean Julian Date of the observations minus 2450000, (p the mean pulsation phase for the observations, 6uD is 
the estimated uniform-disk angular diameter, b the mean baseline of the observations, V^^p the expected value of V'^ for the uniform-disk 
angular diameter, baseline and effective wavelength (696.2 nm), V^^^^ is the mean value of the N observed values of V^, and (V^^p-V^^^^) 
is the difference between the expected and observed values of . Further details are given in the text. 
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for L Car v. s Car. Linear regression fits of the form y = bx 
were made to each of the three plots. 

In the case of s Car v. /3 Car, b — 1.001 ± 0.006 indicat- 
ing that these two cahbrators are consistent with each other. 
However, in the case of t Car v. /3 Car, 6 — 0.898±0.007, and 
for L Car v. s Car, b = 0.897 ± 0.009. Since both these latter 
two plots show essentially the same slope, which differs sig- 
nificantly from the value of unity expected if the calibrators 
were consistent, it suggests that the value adopted for the 
uniform-disk angular diameter of t Car is too small. One 
possibility for this inconsistency is that t Car is a binary 
system but there is no evidence in the literature to support 
this. A careful examination of the SUSI data for the signa- 
ture variations expected from a binary system has proved 
negative. It appears that the value adopted for the uniform- 
disk angular diameter is too small. Although it is hard to 
justify, arbitrarily increasing the angular diameter of t Car 
to 1.77 ± 0.12 mas, ~ 2a greater than predicted, gives b val- 
ues closest to unity for t Car v. (3 Car and t Car v. s Car 
(6 = 1.001 ± 0.011 for L Car v. s Car and b = 1.003 ± 0.009 
for t Car v. /3 Car). The implications of making this change 
to the value of the uniform-disk angular diameter of t Car 
will be discussed in Section |8] but the analysis will initially 
be completed with the data in Table 3] which is based on 
the uniform-disk angular diameter of t Car being equal to 
1.55 ±0.12 mas. 



5.4 Limb-Darkening Factors 

In order to determine the true, limb-darkened angular di- 
ameters, the uniform-disk angular diameters have to be 
multiplied by limb-darkening factors (px) that are depen- 
dent on the centre-to-limb (CTL) intensity distributions for 
the star. The CTL intensity distributions are dependent on 
wavelength, and on the effective temperature (Tefi), surface 
gravity (log (j) and composition ([Fe/H]) of the stellar atmo- 
sphere. 

The form of the CTL intensity distri bution and its 
phas e dependence has been questioned. iMarengo et al.l 
l|2002 . 2003) have computed centre-to-limb intensity distri- 
butions for the Cepheid Gem and shown that their distri- 
butions computed for hydrodynamic models in spherical ge- 
ometry differ from those for hydrostatic, pla ne-parallel mod- 
els, p articularly at certain ranges in phase. iNardetto et al.l 
((2006b) have used a hydrodynamic model of S Cep to derive 
intensity distributions in the continuum and in four spec- 
tral lines. They found that limb-darkening in the continuum 
revealed a systematic shift in phase of the derived angu- 
lar diameter of 0.02. However, the distance is not affected 



because it is linked to the amplitude of the angular diame- 
ter curve, which is only slightly changed by the shift effect. 
They further claim that considering the time-dependence 
of limb-darkening does not seem to be a priority for the 
IBW method. The validity of this claim is considered in 
Section lSTTl and the use of a phase-dependent limb-darkening 
factor is justified for the wavelength of the observations pre- 
sented here (696 nm). 

In the absence of a CTL intensity distribution computed 
for a hydrodynamic model of £ Car, as d one for C Gem, use 
has been made of the tabulation of px bv lDavis et al.l (|200G| ) 
computed for the extensive grid of centre-to-limb intensity 
variations given by Kurucz (1993a b) for his model atmo- 
spheres. 

A v alue of [Fe/H] = 0.3 ha s been adopted for £ Car 
following [Cavrel de Strobel et al.l l|l997h . iTavloil l| 19991 ) has 
tabulated 47 values of T^g and \ogg as a function of phase 
and these data have been adopted to establish a paga versus 
phase curve. For each phase tab ulated by Taylor a value for 
P696 was interpolated from the iDavis et al.l (|200G| ) tabula- 
tion. A sixth order Fourier series was fitted to the result- 
ing values of paga versus phase (the lowest order to give a 
smooth and accurate representation of the data; the data 
has a standard deviation from the fitted curve of 0.00024). 
The coefficients of the fit were then used to compute val- 
ues for Pegs for the phases of the SUSI uniform-disk angular 
diameter determinations. These values are listed in Table |4] 
and plotted in Figure [21 Figure [5] also includes the values 
determined for the Taylor phases and the curve fitted to 
them. Taylor included estimates of the uncertainties in the 
values of Tch and log g in her tabulation and from these un- 
certainties it is estimated that the values of pme in Table [4] 
are accurate to ±0.001, with the caveat that we are assum - 
ing that the Kurucz models upon which iDavis et al.l (|200G| ) 
based their tabulation of limb-darkening factors are appro- 
priate for £ Car. The mean value of p696 for the fitted curve 
is 1.057. 



5.5 The Limb-Darkened Angular Diameters 

The uniform-disk angular diameters have been converted to 
limb-darkened angular diameters {6lo) in Table [4] using the 
listed limb-darkening correction factors. The uncertainties 
in the uniform-disk and limb-darkened angular diameters 
given in Table [4] do not include the uncertainty in the effec- 
tive wavelength, which is a systematic uncertainty of ±0.3%. 
Its effect on the mean angular diameter and on the distance 
to £ Car will be discussed in Section [8] The uncertainty in 
the limb-darkened angular diameter in Table|4]does however 
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Figure 2. Limb-darkening factor P696 v. pulsation phase for 
I Car. Tlie curve is tlie sixtli order Fourier series fit to tlie filled 
cir cles wh i ch are the values determined for the phases tabulated 
bv lTavloil l|l999l V The open squares are the values for the epochs 
of the SUSI determinations of the uniform-disk angular diame- 
ters, computed from the coefficients of the Fourier fit. Further 
details are given in the text. 



include the uncertainty in the limb-darkening correction fac- 
tor. 



6 THE SPECTROSCOPIC DATA 

In order to determine the distance and mean diameter of 
£ Car the radial displacements of the stellar surface as a 
function of phase are required for combination with the 
limb-darkened angular diameters. The starting point for es- 
tablishing the radial displacements is the observed radial 
velocity curve as a function of phase. 



6.1 The Radial Velocity Data 

The most extensi ve radial velocity data set is that by 
iTavlor et all (|l997l ) based on 67 spectra from Mount Stromlo 
Observatory (MSO) in Australia and 70 spectra from Mount 
John University Observatory (MJUO) in New Zealand. 
Early measurements of radial velocity induced spectral line 
displacements were generally ma de by line ce ntroid or bisec- 
tion estimates made by eye (e.g. lDawe|[l969h which was the 
last p ublished velocity curve for I Car prior to lTavlor et al.l 



Taylor et al] l| 19971 ) determined radial velocities by av- 



eraging the radial velocities f or 19 metallic lines me asured 
by the line-bisector method (jWallerstein et al]|l992l ) from 
an average of depths 0.7, 0.8 and 0.9 (continuum at 0.0 and 
core at 1.0). The MSO and MJUO sets of data were ini- 
tially phased using the ephemeris of IShobbrookl l|l992l ) but 
iTavlor et al.l (|l997l ) found a small phase shift between the 
earliest MJUO observations and later observations and de- 
rived a small phase correction to bring all the MSO and 
MJUO observations into phase alignment. We have adopted 
their tabulated data that include the phase adjustments 
with three exceptions. The Julian Dates 2449408.65 and 
2450380.57 for MSO data and 2449408.66 for MJUO data 



correspond to times during daylight hours so data for these 
dates have been omitted leaving a total of 134 data points. 

iBersierl (|2002l ) lists HJDs for 19 radial velocity measure- 
ments for I Car but does not give phases. For I Car the max- 
imum difference between HJD and JD is negligible at less 
than 3 x 10~^ of the pulsation period and has been ignored. 
Phases have been computed with Shobbrook's ephemeris 
and the re sulting data are in good phase agreement with 
the data of ITavlor et al.l l|l997l ) , discussed abo ve, but there 
is an offset in radial velocity. iKervella et al.l (l2004b) first 
noted this and chose 'to shift the Taylor et al.l ( 19971 ) data 



set by —1.5 km s to bring all the data o n the well estab- 
lished CORAVEL s y stem of iBersieil (I2OO2I )'. Since the data 
set of iTaylore t al (1997') is significantly more numerous 
than that of Isersier (2003 ) (134 data p oints versus 19), and 
the decision has been made to adopt the lTaylor et al.l (|l997l ) 
as the basic data set, we have adjusted the Bersier radial 
veloci ties. This decision is supported by the work of iKissI 
(|l998h who found that, while 'CORAVEL measurements 
have excellent internal accuracy, their absolute values are 
very uncertain'. The Bersier data, determined by the cross- 
correlation tech nique rath e r tha n the line-bisector method 
employed by Taylor et al.l (|l997l ) , were combined with the 
Taylor et al. data with a range of offsets for the Bersier ra- 
dial velocities. For each offset a sixth order Fourier series 
was fitted over the range in phase from -0.037 to 0.75 (the 
reasons for the choice of the order of the fit and of the range 
in phase will be discussed in Section [6. 2p and the minimum 
value of corresponds to the addition of 2.0 km s~^ to 
the Bersier data. In view of the resulting good agreement 
between the Bersier and Taylor et al. data over the whole 
pulsation cycle including the amplitude of the radial veloc- 
ity variation, the Bersier data have been accepted in spite 
of the different technique used for measuring the line shifts. 

iPetterson et~al] l|2005l ) have measured 34 radial veloci- 
ti es at the MJUO and while these are more recent than those 
of iTaylor et al] (|l997l ). the technique for measuring the ra- 
dial velocities by the line-bisector method is identical. As 
well as giving the JDs of the observation s the authors have 
listed phases based on an ephemeris by IPell l| 19761 ) which 
are significantly offset from the data of Bersier and Taylor 
et al. We have compu ted the phases fro m the given JDs, 
using the ephemeris of IShobbrookl l|l992l ). and find a small 
phase offset from the Taylor et al. data suggesting a glitch 
or period change between the Taylor et al. and the later Pet- 
terson et al. observations. In order to establish the optimum 
phase adjustment the Petterson et al. data were combined 
with the Taylor et al. data with a range of phase offsets. 
For each offset a sixth order Fourier series was fitted over 
the range in phase from -0.037 to 0.75, as in the case of 
the Bersier radial velocities, and the minimum value of 
corresponds to t he su btraction of 0.02 in phase from the 
IPetterson et al.l (|2005h data. The adjusted Petterson et al. 
data are in excellent agreement with the combined Taylor 
et al. and Bersier radial velocity data with the exception 
of the five data points for JD 2451163.0392, 2450683.1543, 
2450683.8039, 2450684.1304 and 2450684.8105 which aU lie 
several times the quoted uncertainties from the combined 
Taylor et al., Bersier, and Petterson et al. plot of radial ve- 
locity against phase. The Petterson et al. data have been 
accepted with the exception of the five discrepant points. 
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Table 5. Details of the Fourier series fits to the selected radial 
velocity data in three ranges in phase. N is the number of data 
points in each range in phase, Np is the order of the Fourier 
series fit, and Phase{M) is the phase at which the separate fits 
were matched. Further details are given in the text. 



Section 


Phase 
Range 


N 


Np 


Phase(M) 


Rising 


-0.037 to 0.75 


137 


6 




Rising to Middle 








0.710 


Middle 


0.701 to 0.896 


42 


6 




Middle to Falling 








0.890 


Falling 


0.862 to 1.031 


38 


6 




Falling to Rising 








0.015 



iNardetto et al.l (|2006ah have published radial velocities 
for £ Car determined by three different methods of measur- 
ing line displacements. These methods differ from those used 
for the data considered so far and result in a range of ampli- 
tudes of the radial velocity curve bracketing the amplitude 
of the curve determined from the Taylor et al., Bersier and 
Petterson et al. curve with none of the three agreeing with 
it. In view of the d isagreement with the other three sources 
the INardetto et aT.i (,2006a,) data have not been included. 



6.2 The Radial Velocity Curve 

Figure [3] shows the assembled radial velocity data versus 
phase together with the fitted curve. Attempts were made 
to fit a Fourier series to the assembled data b ut even with a 
sixteenth order series the fit was poor in parts. Fl^vloij ((199^) 
overcame this difficulty by fitting a sixth order Fourier fit 
to the ascending branch of the curve and a ninth order fit 
to the descending branch. After some experimentation it 
was found that the good fit shown in Figure |3] could be 
obtained by dividing the curve into three ranges in phase 
and making a separate sixth order Fourier series fit to each. 
The ranges of the fits, which overlap, are listed in Table [5] In 
each case the reduced values indicated that the published 
uncertainties in the radial velocities are optimistic. Adopting 
uncertainties of ±0.45 kms~^, ±0.50 kms~^, and ±1.0 kms^^ 
for the rising, middle, and falling sections respectively gave 
reduced values close to unity. The coefficients of the fits 
were employed to compute radial velocities at intervals of 
0.005 in phase over the ranges of the fits and plots were made 
of the fits in the overlapping regions. The phases for the 
transitions from one fit to the next were chosen by inspection 
of these plots and were the phases at which the overlapping 
fits were closest to each other. For the transition between the 
rising and middle sections the radial velocities at phase 0.710 
were adjusted to an intermediate value to give a smooth 
transition between the fits on either side. Similarly, for the 
transition between the middle and falling sections, the radial 
velocity at phase 0.890 was adjusted and, for the transition 
between the falling and rising sections, the radial velocity 
at phase 0.015 was adjusted. The adjustments to the radial 
velocities were small, averaging ~0.3% of the mean radial 
velocities at the transitions. The fits are summarised and 
the phases at which they were matched are listed in Table (5) 
The radial velocity curve shown in Figure |3] assem- 
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Figure 3. The assembled radial velocity curve v. phase for 
£ Car. T he line is the fit to the data; filled circles are MSO 
data fromlTavlor et aL open circles a re MJUO data from 

iTavlor et al.l l|l997h : open squares ar e data byllersieilliooi); and 
crosses are data by iPetterson et al.l ll2005h . Details are given in 
the text. 

bled from the Fourier series fits listed in Table O has been 
adopted for the subsequent analysis 

6.3 The p-Factor 

The radial velocity curve in Figure |3] is for the measured ra- 
dial velocities but does not represent the true radial velocity 
of the stellar surface since it is a value integrated over the 
stellar surface. It includes projection and limb-darkening ef- 
fects that vary from the centre to the limb of the star. These 
effects depend on spectral line shape, which is not only a 
function of phase, but is also affected by the velocity struc- 
ture within the line-forming region and the contributions 
from the different layers of the atmosphere. Conventionally, 
in the absence of sufficiently detailed modelling of the stel- 
lar atmosphere, the measured radial velocities are multiplied 
by a constant, known as the p-factor, to correct them to the 
radial velocity of the stellar surface. There are a number of 
issues to be considered in the choice of p since any error or 
uncertainty in the value adopted translates directly to the 
distance determined for the Cepheid when the spectroscopic 
and interferometric data are combined. 

The evaluation of p-factors is genera lly derive d from 
model stellar atmosphere models but Mcr and et al.l ()2005h 
have presented a measured value of 1.27±0.06 for the p- 
factor for 5 Cep based on interferometric measurements with 
the C HARA Array and on the Hubble Space Telescope par- 
allax (|Benedict et al.|[2002t ). The accuracy is limited by the 
parallax and the authors conclude that theoretical studies 
using realistic hydrodynamical codes are needed. 

The literature dealing with the theoretical evaluation of 
the p- factor from model atmosphere studies is extensive but 
with no clear results applicable to £ Car. While lSabbev et al.l 
(1995) claim that the phase dependence of p increases the 
Baade-Wesselink (BW) radius by ~ 4 — 6% , depending on 
the co nstant value of p used for comparison, INardetto et al.l 
claim that their choice of a constant p-factor for the 
Interferometric Baade-Wesselink method (IBW), compared 
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to a time-dependent one, leads to a systematic error of the 
order of only 0.2% in the final distance determination for 
5 Cep. 

The projection factor is also sensitive to the centre-to- 
limb intensity distribution or limb darkening, but it is not 
clear whether a mean value for the limb-darkening is ade- 
quate or whether a phase dependence is significant for the 
p-factor. However, it is clear that the phase dependence of 
limb darkening is significant for the conversion of uniform- 
disk angular diameters to limb-darkened angular diameters 
at 696 nm, as mentioned in Section 15.41 and justified in Sec- 
tion [H] 

In most cases more detailed modelling and, in particu- 
lar, the hydrodynamic modelling in spherical geometry that 
predicts a phase dependence of p has been done for a partic- 
ular Cepheid, most commonly 5 Cep or C, Gem, and the re- 
sults are not in a form that can be scaled to other Cepheids. 
The concluding recommendation of these studies is generally 
that each Cepheid should be individually modelled. This is 
not within the scope of our programme and we have there- 
fore decided to use a fixed value but to make available all 
the relevant data to enable the results to be updated by oth- 
ers when improved values for p, or an appropriate model for 
I Car have been developed. 

A brief summary and discussion of the fixed val- 
ues of p adopted for I Car in t he literature is appr opri- 
ate at this point. As noted by iTavlor et al.l (|l997f ) the 
value of p depends not only on the technique used to 
measure the radial velocities, but also upon the strength 
of the lines chosen and their wavelengths. Ta ylor ot al. 
lll99it adopted a value of 1.3 8±0.03 ifAlbrow fc Cottrelll 
1 19941 ') but revised their results jTavlor fc Boot h' '1998^ us- 
ing p = 1.39 - 0.03 log P = 1.34 dHindslev fc Bell 1986; 
iGieren. Barnes fc Moffettlll993l ). fKervella et al.! (|2004bl in 
a comparison of the interferometric and surface brightness 
techniques, used the same formula with p = 1.343 for £ Car. 
In this comparison they only considered ra dial velocity 
points in the phas e interval 0.0 to 0.8 foll owing IStorm et al.l 
(|2004l ). However, iKervella et all (j2004al ) adopted a value 
of 1.36 for p for all the 7 Cepheids in their programme 
which included t Car and j ustified it on the grounds that 
iBurki. Mavor fc Bend (|l982l 'l had shown that this v alue was 
appro priate for the radial velocity measurements bv lBersieil 
(|2002 >') that they had used. iNardetto et al.l (|2007l ) have used 
a hydrodynamic model of I Car to validate a spectroscopic 
method of determining the p-factor in which it was divided 
into three sub-concepts. While their work is not directly ap- 
plicable to our data, because it was restricted to a single 
specific spectral line and the method employed for measur- 
ing the line differed from that for th e data we are using, 
they derived a value of 1.27±0.02 for p. ICroenewegenl (|2007l ) 
has evaluated a relationship between the p-factor and pulsa- 
tion period based on five Cepheids with interferometrically 
measured angular diameter variations and known distances 
taken from the literature. Based on a total of seven stars 
with periods in the range 5-35 days it is claimed that there is 
no evidence for a period dependence of the p-factor although 
values found range from 1.193 to 1.706, albeit with large un- 
certainties in the majority of cases. For I Car the p-factor is 
found to be 1.193±0.058 ±0.120 based on the distance of 
498±50 pc determined by [Benedict et all (120071 ) with the 
Hubble Space Telescope Fine Guidance Sensors. The first 



uncertainty is from the fitting process and the second is due 
to the uncertainty in the distance. The majority of distance 
determinations to I Car to date, which are listed in Table [S] 
and include the distance determined here, suggest a larger 
value for the distance to I Car, implying a larger value for the 
p-factor of the order of 1.29 if determined by the approach 
employed by Groenewegen. Based on the five Cepheids anal- 
ysed, which had interferometrically determined angular di- 
ameters and distances, Groenewegen (2007) concluded that 
a constant value of p = 1.25 ± 0.05 was appropriate. But, 
as has been discussed, at least in the case of I Car a larger 
value is indicated. 

To summarise, there is no agreement on the optimum 
fixed value of the p-factor to use although it is clear that 
it depends on the individual Cepheid and on the details of 
the observing and analysis techniques used for both spec- 
troscopy and interferometry. Faced with these difficulties we 
have adopted a constant value of 1.30±0.05. The influence 
of this decision on the mean diameter and distance of I Car 
will be discussed in Section [S] 



6.4 The Radial Displacement Curve 

Integration of the radial velocity curve shown in Figure [3] 
gave the radial velocity of the centre of mass V-y equal 
to -|-4.13±0.01 kms~^ where the uncertainty has been es- 
timated using the bootstrapping metho d. This differ s from 
the value of -|-4.21±0.01 kms~ given bv lTavlor etHI (|l997D 
but is based on the assembly of a larger body of radial ve- 
locity data. 

The radial displacement of the stellar surface in solar 
radii AR{(j)), as a function of phase </!>, was found by inte- 
grating the radial velocity Vrv (<?!>), after correction for V^, 
using 



Ai?(0) 



pP 



(Frv(</>) - n)d0 



(4) 



where P is the pulsation period in seconds and Rq is the 
solar radius in km. 

The integrations were made at intervals in phase of 
0.005 and the resulting values of radial displacement, af- 
ter the determination and subtraction of the mean value, 
were fitted with a Fourier series. Series of increasing order 
were fitted and it was found that the reduced value for 
the fits had a minimum value for a sixteenth order fit. This 
has been adopted for the subsequent analysis. The uncer- 
tainty in the radial displa cements has been evalua ted fol- 
lowing iTavlor et al.1 (Il997l ) using the expression by iBalonal 
1 1977 



o-(RD) = 



pPa(RV) 



(5) 



where a(RD) is the uncertainty in radial displacement in so- 
lar radii, P is the pulsation period in seconds, a(RV) is the 
standard deviation of the radial velocities about the fitted 
curve (ct(RV) = 0.6 km s~^), Rq is the solar radius in km, 
and N is the number of observations. Substitution in equa- 
tion ((U gives cr(RD) = 0.13 Rq (less than 0.4% of the total 
radial displacement due to the Cepheid pulsation). 
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7 THE COMBINATION OF 
INTERFEROMETRIC AND 
SPECTROSCOPIC DATA 

The relationship between the limb-darkened angular diame- 
ter ^LDobs and the radial displacement of the Cepheid surface 
is given by 



+ 9.298 



\ d 



(6) 



where 0ld is the mean limb-darkened angular diameter 
(the limb-darkened angular diameter at zero displacement), 
AR{4)i) is the radial displacement in solar radii for the ith 
observation at phase 0^, d is the distance in pc, and the con- 
stant converts the term in parentheses into ma s using valu es 
for the solar radius and parsec in metres from lCoxl l| 19991 ). 

A weighted linear least squares fit has been made to 
OLoi(f>i) versus AR{(j)i) to determine Sld and d for ^Car. 
A small phase shift was found between the angular diam- 
eter and radial velocity data, which is not surprising since 
the epoch of the first SUSI observation was more than 40 
pulsation periods after the last radial velocity observation. 
The phase offset was established by repeating the fit with a 
range of phase offsets to find the minimum value of reduced 
for the fit. The phase offset was found to be a correction 
of -0.0635±0.015 to the phases of the radial displacements 
relative to the phases of the angular diameters. The fit using 
this offset gives the mean limb-darkened angular diameter 
of £ Car equal to 2.981±0.005 mas and the distance to £ Car 
equal to 523±15pc. The observational data are shown in 
Figure |4] with the fitted line. 

In Section 15.31 the apparent inconsistency of the 
uniform-disk angular diameter of t Car was discussed and 
the possibility of increasing it to 1.77±0.12mas was pro- 
posed although all the analysis to this point has been 
carried out with our estimate for its angular diameter of 
1.55±0.12 mas. In order to investigate the effect of in- 
creasing the angular diameter the analysis to this point 
has been repeated with the angular diameter of t Car 
equal to 1.77±0.12 mas. The results are: phase offset - 
0.064±0.001, mean limb-darkened angular diameter of £ Car 
equal to 2.999±0.005 mas and the distance to £ Car equal to 
526±15pc. 

At this point it is stressed that the uncertainties are 
those given by the fits. However, it can be seen that the 
change in the angular diameter of t Car has essentially neg- 
ligible effect on the distance determination compared with 
the uncertainty in the distance, but does increase the mean 
angular diameter by more than three times the formal un- 
certainty from the fits. The difference in the mean angular 
diameter and systematic uncertainties that have not been 
considered so far are discussed in Section [S] 

The limb-darkened angular diameters versus phase for 
£ Car, assuming the angular diameter of t Car equal to 
1.55±0.12 mas, are shown in Figure[5]together with the limb- 
darkened angular diameter versus phase curve computed for 
the fit of limb-darkened angular diameter versus radial dis- 
placement. 
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Figure 4. The limb-darkened angular diameter versus radial dis- 
placement for £ Car. The filled circles with error bars are the 
observed values of limb-darkened angular diameter from Table |4] 
and the line is the fit to angular diameter versus radial displace- 
ment using equation ([6J. Further details are given in the text. 
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Figure 5. The limb-darkened angular diameter versus phase for 
£ Car. The filled circles with error bars are the observed values 
of limb-darkened angular diameter from Table |4] and the line is 
the computed limb-darkened angular diameter versus phase from 
the fit to angular diameter versus radial displacement using equa- 
tion (|6]l. The data points and curve in the phase interval 0-0.5 
have been repeated for the phase interval 1.0—1.5. Further details 
are given in the text. 



8 DISCUSSION 

In Section [7| two values were presented for the mean limb- 
darkened angular diameter of £ Car: 2.981±0.005 mas for the 
uniform-disk angular diameter of the calibrator l Car equal 
to 1.55 mas and 2.999±0.005 mas for t Car equal to 1.77 mas. 
Increasing the uniform-disk angular diameter of t Car pro- 
duced greater consistency between the three calibrators, as 
discussed in Section [5.31 but the change was arbitrary and 
we cannot be certain that the larger value is correct. We 
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therefore adopt the mean of the two results and assign an 
uncertainty to cover the uncertainties of the two values. In 
addition there is a systematic uncertainty due to the un- 
certainty in the effective wavelength of ±0.3 per cent. Thus 
the final value for the mean limb-darkened angular diameter 
of I Car, with the systematic uncertainty in parentheses, is 

= 2.990±0.014 (±0.009) mas. This value is in excellent 
agreement with the only ot her direct interferometric deter- 
mination bv iKervella et al.1 (t2004a) of 2.988±0.012 mas. 

The corresponding two values for the distance presented 
in Section[7l namely 523±15pc and 526±15pc, are in good 
agreement and we adopt 525±16pc for the distance. There 
are two systematic uncertainties in the distance due to 
the effective wavelength (±0.3 per cent) and the p-factor 
(±3.8 per cent). The value for the distance to I Car, with the 
systematic uncertainty in parentheses, is 525±16 (±20) pc. 
Combining the uncertainties quadratically gives the distance 
as 525±26 pc. 



8.1 A Constant versus a Phase-Dependent 
Limb-Darkening Factor 

The limb-darkening factor for converting uniform-disk an- 
gular diameters to limb-darkened angular diameters is ex- 
pected to vary with pulsation phase and we have taken this 
into account as discussed in Section [5.41 However, the ques- 
tion has been asked as to whether this was necessary and 
would a constant value, equal to the mean value taken from 
the curve in Figure [2j have given different results. To ex- 
amine this, the entire analysis has been repeated using the 
mean value of the limb-darkening factor of 1.057, given in 
Section [5.41 in place of the phase-dependent values given in 
Tabled 

The mean limb-darkened angular diameter was found to 
be unchanged as expected. The mean limb-darkened angular 
diameter of I Car was 2.979±0.017 mas for the uniform-disk 
angular diameter of the calibrator l Car equal to 1.55 mas 
and 2.997±0.017 mas for t Car equal to 1.77 mas, where the 
statistical and systematic uncertainties have been combined 
quadratically. In each case the values are slightly smaller 
than the values for the phase-dependent limb-darkening fac- 
tor (0.002 mas) but this is a result of the rounding of the 
mean value of the limb-darkening factor. 

The distance to I Car is changed because the mean limb- 
darkening factor results in a smaller angular pulsation am- 
plitude (0.525±0.018 mas compared with 0.560±0.018 mas). 
The corresponding values for the distance are 559±17pc and 
561±17pc compared with 523±15pc and 526±15pc for the 
phase-dependent limb-darkening factor. The mean of the 
two values in each case is 560±18pc and 525±16pc. The 
difference is significant and just ifies the use of the ph ase- 
dependent limb-darkening factor. IKervella et al.l (|2004bl ) ne- 
glected the phase dependence of the limb-darkening factor 
in their s tudy of I Car. Their d ecision was based on the esti- 
mate bv lMarengo et al] (|2003| ) that the variation would be 
less than 0.3 per cent peak to peak in the H band for Gem 
and the fact that it would be even less in the K band. At 
696 nm the limb-darkening factor varies by more than 1.3 per 
cent, as shown in Figure |2l and cannot be ignored. 



8.2 A Comparison of Distance and Radius 
Determinations 

Table [6] lists recent values for the distance to I Car from the 
literature together with the value determined in this work for 
a phase-dependent limb-darkening factor. For the latter the 
statistical and systematic un certainties have b e en com bined 
quadratically. The value by iTavlor fc BoothI (Il998l) using 
the B E method, which succeeds the value by iTavlor et al.l 
([1993), has an unrealistic published uncertainty of ±4 pc be- 
cause systematic effects have no t been taken into account. 
ICieren. Fouque fc Gomej (|l997l ) have shown that a system- 
atic uncertainty of the order of ±3 per cent should be taken 
into account for the surface brightn ess method and the un- 
certainty for the lTavlor et al.l (|l997l ) distance in Table[6]has 
been increased to include this additiona l uncert ainty. It is 
noted that the values by IKervella et al.l (l2004bl) using the 
IBW and IRSB methods! and bv iBarnes et all (|2005l ) by 



the IRSB method, do not include an allowance for the sys- 
tematic uncertainty in the p-factor adopted. 

The distance values in Table |S] that involve the p-factor 
have used different values to that adopted here and, while it 
is recognised that different values may be appropriate for dif- 
ferent methods, a column has been added to the table with 
the distances scaled to a common value of p = 1.30±0.05. 
The uncertainties associated with this column have been 
increased to combine quadratically the uncertainty in the 
p-factor with the published uncertainty where appropriate. 
The uncertainty in t he value by Taylor e t al. (1997) has been 
increased further as lGieren et al.l (|l997l ) included a smaller 
uncertainty in the p-factor (±2.5%) in their evaluation of the 
systematic uncertainty. Examination of the distances for the 
common p-factor show excellent agreement for all values. 

Columns for the mean radius of I Car are included in 
Table |6l one for the published values and one for a common 
p-factor, following the same approach as for the distance. 
The uncertainty in the limb-darkened angular diameter for 
the IBW determinations has been taken into account al- 
though it is negligible compared with the uncertainty in the 
distance. Examination of the mean radii for the common 
p-factor also show excellent agreement for all values. 

The mean total range in limb-darkened angular diame- 
ter has been taken from the curves fitted to limb-darkened 
angular diameter versus radial displacement for the two val- 
ues of uniform-disk angular diameter for the calibrator l Car. 
It is 0.560±0.018 mas, which corresponds to 18.7±0.6 per 
cent of the mean stellar diameter. 



9 CONCLUSION 

The measurement of the angular pulsations of Cepheid vari- 
ables and the combination with spectroscopically deter- 
mined radial displacements of the stellar atmospheres was 
one of the key programmes for which SUSI was developed. 
Only recently has adequate sensitivity been achieved for this 
programme to be undertaken and the first results are those 
for I Car presented here. 

The angular pulsation of I Car has been measured 
with good cover in pulsation phase and careful attention 
has been paid to uncertainties in the measurements, in 
the adopted calibrator angular diameters, in the projected 
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Table 6. The distance to £ Car. p is the p-factor, the column headed ID contains acronyms for the methods employed with the key at 
the foot of the table, and the key to the numbered references is also at the foot of the table. Further details are given in the text. 



Distance Radius p ID Distance Radius Reference 

for p = 1.30 for p = 1.30 
(pc) (Rq) (pc) (Rq) 



550±17 


173±5 


1.34 


BE 


534±22 


168±7 


1 




182^^ 


1.343 


IBW 




17611° 


2 


560±23 


179±7 


1.343 


IRSB 


542±30 


173±10 


2 


559±19 


179.9±6.4 


1.343 


IRSB-Bay 


541±27 


174±9 


3 


485±64 






RH 






4 


498±50 






HST 






5 


525±26 


168.8±8.2 


1.30 


IBW 


525±26 


169±8 


This work 



Key to ID acronyms: BE - Barnes-Evans; IBW - Interfcromctric Baadc-Wcsselink; IRSB - Infrared Surface Brightness; IRSB-Bay - 
Infrared Surface Brightn ess ( Baycsian); R H - Re vised Hipparcos Parallax; HST - Hubble Space Tele scop e Fine Guidance Sen sors. 
Key to References: 1 - [Tavlor fc BoothI lll998l '); 2 -|Kcrvclla ct al. (2004M; 3 - iBarnes etahl ||2005|) ; 4 - Ivan LeeuwenI teOOTD ; 5 - 

iBenedict et al.. (.200^ . 



value of at zero baseline, and to systematic eflFects. A 
phase-dependent limb-darkening factor, to convert uniform- 
disk angular diameters to limb-darkened angular diame- 
ters, was established based on Kurucz model atmospheres 
and its use, rather than a fixed mean value for the limb- 
darkening factor, was justified for 696 nm. The resulting 
mean limb-darkened angular diameter is 2.990±0.017 mas 
(i.e. ±0.6 per cent) with a maximum-to-minimum amplitude 
of 0.560±0.018 mas corresponding to 18.7±0.6 per cent in 
the mean stellar diameter. The projected value of at zero 
baseline shows no evidence at 696 nm of th e circumstellar 
envelo pe observed in the A'^ and K bands bv lKervella et all 
l|2006t ). 

A radial displacement curve has been computed from 
carefully selected radial velocity measurements from the lit- 
erature and this has been combined with the interferometric 
data to determine the mean radius and distance to i Car. 
The value for the distance is 525±26 pc and the mean radius 
is 169±8iZQ. These values have been compared with values 
in the literature and excellent agreement is found, particu- 
larly when all values are reduced to a common p-factor. 

The uncertainty in the p-factor is a major contributor 
to the uncertainties in the distance and mean radius. The 
SUSI data may be useful in the future if £ Car specific hydro- 
dynamic spherical models are generated that consider the 
possibility of a phase-dependence of li mb-darkening and the 
p-fact or as discussed, for examp le, bv lSasselov fc Karo vskal 
1 19941 ) and lMarengo et al.1 l|2002l '). To facilitate such possible 
applications the calibrated values may be obtained from 
the lead author (JD). 
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